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Abstract ■ The analysis of sea level time 
variations recorded by shore tide gauges for 
evidence of recent and contemporary vertical 
crustal movement is an established technique. 
However, the additional refinement of first 
modelling the numerous physical influences on 
sea level — which are of considerable importance, 
but additional complexity in the coastal zone — 
has not been fully exploited. 

In this paper, attempts to specify and eval- 
uate such physical influences are reviewed with 
the intention of identifying problem areas and 
promising approaches. For a routinely viable 
procedure, it is considered important to limit 
the input data, which would be needed on a cont- 
inuing basis to model physical phenomena, to 
widely and regularly available data records. 

An example of linear modelling based on air/ 
water temperatures, atmospheric pressure, river 
discharges, geos trophic and/or local wind veloc- 
ities, and including forced period terms to allow 
for the long period tides and Chandlerian polar 
motion is evaluated and applied to monthly mean 
sea levels recorded in Atlantic Canada. Refine- 
ment of the model to admit phase lag in the 
response to some of the driving phenomena is 
demonstrated. Spectral analysis of the residuals 
is employed to assess the model performance. The 
results and associated statistical parameters are 
discussed with emphasis on ellucidating the sens- 
itivity of the technique for detection of local 
episodic and secular vertical crustal movements, 
the problem areas most critical to the type of 
approach, and possible further developments. 

Introduction 

Vanicek [1978] has described a method whereby 
some of the known physical influences on sea 
level can be removed by a process of mathematical 
modelling. He concludes that, with such a model 
restricted to the effects of atmospheric pressure 
and temperature, river discharge, and three basic 
tidal constituents: "...we should be able to 
detect possible local episodic vertical movements 
(of duration more than A months) of magnitude of 
10cm and more...". It is further stated that, 
by extending and refining the model to include 
additional known effects, it should be possible 
to see movements down to 5cm in the sea level 
record. 

The work discussed in this paper is merely a 
continuation and extension of Vanfceks approach. 
In particular, the basic model has been substant- 
ially retained, but refined by introducing an 
air-to-water heat transfer model with provision 
for time lag, and also by including wind stress 
terms based on either the observed wind vector 
or the geos trophic wind vector derived from the 
observed atmospheric pressure field. 

So far, only results of preliminary and test 
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computations are available, based on data from 
eight locations in the Canadian Maritime provinc- 
es (fig. 1). However, these results have been 
adequate to indicate some strengths and weakness- 
es in the basic approach and to suggest areas 
where immediate refinement could be advantageous. 
Also, some model improvements with the best 
potential for future development have been high- 
lighted. 

The Problem and Proposed Approach 

A detailed description of the problem, the 
proposed approach, and the reasons for choosing 
such an approach have already been set out by 
Vanfcek [1978]. It will suffice here to provide 
a brief summary. 

The problem is to devise a viable technique 
which would enable use of sea level records, 
derived from tide gauge observations, to detect 
contemporary vertical crustal movements over 
periods from a few months to several years. Both 



Figure 1. Map of Canadian Maritime provinces 
showing location of tide gauges 
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secular movements (i.e. linear trend) and irreg- 
ular or episodic movements are sought. To achie- 
ve sufficient accuracy in isolating these phenom- 
ena, a variety of other known physical influences 
on the sea level record must be removed. It is 
proposed that this be done by empirically model- 
ling the known effects, so that much of the 
variance in the sea level record on that account 
can be eliminated. 

Two critical factors are inherent in this 
approach. Firstly, the tide gauge environment 
may be substantially more complex than that which 
is experienced in the deep oceans. For example, 
esturine topography, the dynamics of coastal 
bathymetry, seasonal freezing of coastal waters, 
and dynamic river discharge influences are some 
of the many factors which may interact with the 
meteorological and tidal effects to anomalously 
perturb the coastal sea level record. Consequent- 
ly, many of the well proven modelling techniques 
devised for the deep ocean environment and much 
of the data observed there may not necessarily be 
extrapolated with any validity to the tide gauge 
environment. Secondly, the method of modelling, 
to be viable for routine monitoring of vertical 
movements, should rely only on readily available 
long-term data records, collected on a regular 
basis. This is a significant and frustrating 
limitation. 


The Basic Model and Method 

Basically, 'the model adopted (table 1) is a 
linear combination of the physical effects to be 
represented, so that the method might be referred 
to as a multiple linear regression analysis. 
However, certain components of the basic model 
may be represented by time series data which has 
been derived from additonal observations through 
sub-models. This is specifically so in the case 
of the water temperature and wind stress terms. 
The sub-models are not necessarily linear. 

A solution for the unknown parameters, along 
with estimates of their variance and a correla- 
tion matrix, is obtained by the conventional 
least squares method. Actual implementation of 
this solution, however, is just a by-product of 
the optimum least squares spectral analysis 
algorithm [Wells and Vanfcek, 1977] which is 
employed. This algorithm, based on the technique 
developed by Vanfcek [1971], primarily provides 
the spectrum of the residual time series, R(t i ) . 
Two particular advantages claimed for this algor- 
ithm are: (a) suppression of the regression model 
components produces no movement of the inherent 
spectral peaks; and (b) equally spaced time 
series data is not mandatory. Spectral analysis 
of the residuals provides a supplementary quant- 
itative assessment of the success of the regress- 


The Linear Regression Model 

S(ti) = C A + C^t* + CpdPfq) + C T 6T(t t ) + C^D^) + C Wt 6W t (t i ) + 


+ t A, C.Oi (u.tj - 6.) + R(tjj) 
j=l J J 


Unknown Parameters 


C A datum bias 


Observed Time Series 


(functions of time t^) 


Cl linear trend coefficient 

Cp air pressure coefficient 

Cq- surface water temperature coeff- 
icient 

Cp river discharge coefficient 

C W( . shoreline tangential component of 
wind stress coefficient 

Cy shoreline normal component of wind 
stress coefficient 

Aj amplitude of five periodic tidal 
components 

0j phase angle of five periodic tidal 
components 

R sea level residual time series 


S monthly mean sea level at tide gauge 

6P monthly mean air pressure variation 
(with respect to temporal mean) 

6x monthly mean surface water temperature 
variation 

5D monthly mean (combined) river dis- 
charge 

6W t shoreline tangential component of 
observed or geostrophic wind stress 

6W n shoreline normal component of observed 
or geostrophic wind stress 

ii>j assumed constant frequency of five 

periodic tidal components as follows: 
u)L = 6 mths. (semi-annual), ii>2 = 12 
mths. (annual), 013 = 14.33 mths. 
(chandlerian or polar), = 8.847 yrs. 
(lunar perigee), 0)5 = 18.613 yrs. 

(lunar nodal). 


Table 1. The linear regression model 




ion model and provides diagnostic evidence by 
depicting the distribution of the remaining var- 
iance. 

Refinement of the Model 

Assuming the desirability of a linear model, 
two further possibilities remain available for 
improvement. The sub-models may be reformulated 
to better portray the physical processes or the 
total linear model may be extended to include 
additional terms, representing known or suspected 
physical effects. Both methods were employed, as 
exemplified respectively by the empirical sub- 
modelling of water temperature and the addition 
of wind stress terms to the basic model. 

A Water Surface Temperature Model . While ' 
thermal variations have been cited as an import- 
ant contribution to sea level changes, earlier 
analyses of the Maritimes data, using a simple 
linear term based on air temperature, produced 
unexpectedly small and weakly determined regress- 
ion coefficients. It was suggested that this 
could be the result of two factors associated 
with the use of air, rather than water, temperat- 
ures. Clearly, the primitive model ignored time 
lag between air temperature variations and their 
consequent effect on water temperatures and thus 
the sea level. Also, air temperatures usually 
fall below freezing point during Maritime winters 
but not all of the gauge locations actually 
freeze, which suggests that a better thermal 
model, capable of accounting for the "winter non- 
linearity" should be constructed from water temp- 
eratures. Unfortunately, regular long-term 
records of water temperature are not readily 
available at most tide gauges. However, the 
viability of the method could be preserved by 
constructing an empirical sub-model for the air/ 
water heat transfer process, based on long-term 
air and water surface temperature records avail- 
able for the Halifax, N.S. tide gauge. The 
assumption might then be made, and possibly 
tested, that such a sub-model could be validly 
extrapolated to the other gauges. Accordingly, 
a linear heat transfer sub-model, incorporating 
time lag, was postulated as follows: 

T(t ± ) = hT(t ± + fit) + k, (1) 

where t(t^) is the predicted water surface temp- 
erature at time t^, T(tg + fit) is the observed 
air temperature at time t^ + fit, fit is a time 
lag, h is a heat transfer coefficient, and k is 
a constant temperature bias. A least squares 
solution for the parameters fit, h, and k was 
obtained using the monthly mean air and water 
surface temperature series available for Halifax. 
Since the observed time series comprise discrete 
values, the admission of a time lag necessitates 
interpolation of the air temperature data. Though 
the simplicity of polynomial interpolation was 
attractive it was found, by experiment, to be 
considerably less reliable than fourier interp- 
olation; a not entirely unexpected outcome when 
the "sine wave" form of the air temperatures is 
considered. When a fourier transformation is 
introduced into the model (1), the observation 


equations become, after linearization: 

_ m 

T(t-j+fit)Ah - [h J jt)AjfilJl( jo)(t^+5t) - 0j^]At 

+ k + [hTCti+fit) - T( ti )] = v ± (i=l,n); (2) 

where fit and h are estimates of the parameters 
fit and h respectively, so that fit = fit + At and 
h = fi +-Ah; Aj and 8j are respectively amplitude 
and phase fourier coefficients of the air temp- 
erature time series; id is the harmonic frequency; 
v^ are the discrepancies in the observations; m 
is the number of harmonics; and n is the number 
of observations. 

Least squares estimates of the parameters and 
their standard deviations are listed in table 2, 
along -with other- relevant information. The 
estimates are apparently well determined, mainly 
because of the large number of degrees of freedom, 
and the overall RMS error of prediction — at about 
1°C — is quite satisfactory. Distribution of the 
residuals was significantly normal. Use of these 
values in the sub-model (1) also requires fourier 
interpolation of the air temperature series. 
Extrapolation to other locations has not yet been 
attempted, as additional water surface temperat- 
ure data is being sought so that the spatial 
coherence of the coefficients fit, h, and k can be 
tested. 

Wind Stress . Two methods of modelling the 
wind stress effect on sea level have been consid- 
ered. The first relies on a sub-model for the 
geostrophic wind component, derived from the air 
pressure gradient field as follows: 
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where Wg and W^ are the east and north components 
respectively of the geostrophic wind; R is the 
earth's radius; p is air density; to is the 


Coefficient 

Least Squares 
Estimate 

Standard 

Deviation 

fit (mths.) 

-0.66 

0.03 

h 

0.628 

0.009 

k (°c) 

3.05 

0.12 

Data 



Span: Jan. 1927 to Nov. 1973 

excluding Jul. 

1933 to Jan 1946 
Number of observations, n = 412 

Number of fourier coefficients, m = 206 

Statistical analysis 

RMS residual = 1.035°C 
Maximum residual = 4.4°C 
Variance ratio = 1.06 



Table 2. Air /water heat transfer model 
coefficients for Halifax, N.S. 
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earth's angular velocity of rotation; and 3P/3A, 
3P/3<£ are the components of the barometric press- 
ure (P) gradient with respect to longitude (A) 
and latitude ( <f>) . Even though this model may 
only account for part of the actual wind field, 
it has been introduced for use at gauge locations 
where suitable observed wind records are not 
available. The suitability of the wind records 
is stressed here because, in many instances, 
monthly mean magnitudes and directions are seper- 
ately compiled from the hourly mileages, without 
proper vector summations. Such data is thus not 
suitable for modelling wind stress effects. How- 
ever, early attempts using the geostrophic wind 
in the Maritimes region — resolved into components 
tangential and normal to the local shoreline and 
squared — have not been particularly successful 
[Vanfcek, 1977]. It is now suspected that these 
results may have suffered from computational 
difficulties and the present author is attempting 
to rectify the analyses. 

As a means of verifying the significance of 
the wind stress effect, the second method — which 
relies on the observed monthly mean wind vector — 
has been tested at Halifax. Although the Halifax 
wind records are burdened with the unsuitable 
compilation techniques previously mentioned, it 
was possible to recompute, without undue effort, 
correct monthly vector averages since the records 
included monthly mean wind mileages within a 
number of "direction sectors". Prior to 1967 the 
horizon was divided into 8 such sectors of 45° 
each, but since that date 16 sectors have been 
employed. North and east components (Wjj and Wg) 


of the observed monthly mean wind vector were 
thus, easily derived. These components are 
employed in the wind stress sub-model thus: 

(Wjj coA a + We i-Ln a ) ^ 

(W£ C.Oia- Wn 4-in a) | (W E C04 a - Wfj 4-tn a) | 

(4) 

where W t and W n are, respectively, tangential and 
normal to shoreline components of the wind stress 
and a is the shoreline azimuth. 

Results 

Results for four of the eight gauge locations 
are displayed in table 3 and figures 2 and 3. 

The analyses of the remaining gauge records (i.e. 
Point Sapin, Pictou, Boutilier Point, and Yar- 
mouth) are not included here because the use- 
fulness of the results is considerable diminished 
by the short data spans — all less than 9 years. 

In some cases, these additional results provide 
evidence of spatial coherence of the phenomena 
studied. This, as well as any significant dis- 
agreements with the results presented, will be 
mentioned in the relevant discussion below. 

Sea Level Data 

The first four lines of table 3 quantitatively 
describe the observed sea level time series S(ti) 
for the selected tide gauges and they are depict- 
ed graphically in figure 2. The disjoint nature 
of the Pointe au Phre record should be noted. 


Location of Tide Gauge 

Pointe au P&re 

Charlottetown 

St. John 

Halifax 

Number of observations 

422 

370 

480 

564 

Number of years 

47 

31 

47 

47 

Epoch 

1927-1974 

1943-1974 

1927-1974 

1927-1974 

Number of gaps in data 

7 

2 

5 

0 

Datum bias, C A (era) 

27.6 * 0.9 

158.6 1 0.4 

426.4 + 0^5 

110.1 ± 0.3 

Linear trend Cl (cm/cent.) 

7.6 * 2.3 

27.0 1 2.2 

37.6 i 1.6 

43.4 t 1.1 

Tide: Amplitude Aj (era) 





Semi-annual 

2.2 10.4 

0.9 1 0.3 

2.1 i 0.3 

0.9 + 0.2 

Annual 

1.5 t 0.5 

4.6 i 0.4 

1.4 i 0.7 

4.0 ± 0.3 

Chandlerian 

0.7 ? 0.4 

0.7 ! 0.3 

1.3 i 0.3 

0.6 1 0.2 

Lunar perigee 

1.6 t 0.4 

0.7 * 0.3 

0.3 i 0.3 

0.2 ± 0.2 

Lunar nodal 

0.5- * 0.4 

1.6 1 0.3 

1.2 1 0.3 

0.3 ± 0.2 

Tide: Ascending node date, Nj (yrs) 





Semi-annual 

1926.81 ± 0.02 

1926.83 ± 0.02 

1926.81 i 0.01 

1926.81 ± 0.02 

Annual 

1926.39 * 0.05 

1926.64 t 0.04 

1926.80 i 0.04 

1926.69 i 0.02 

Chandlerian 

1926.66 .♦ 0.10 

1926.49 t 0.07 

1926.56 J. 0.04 

1926.52 i 0.06 

Lunar perigee 

1932.73 i 0.38 

1932.80 ± 0.55 

1930.28 i 1.13 

1926.59 ± 1.39 

Lunar ’ nodal 

1936.50 > 3.00 

1930.16 ± 0.52 

1933.54 + 0.70 

1928.21 ± 2.31 

Temperature, C T (cm/°C) 

-0.15 ± 0.09 

-0.033 1 0.087 

0.22 i 0.07 

0.004 i 0.076 

Air pressure, Cp (cm/rabar) 

-0.74 ± 0.11 

-0.90 i 0.06 

-0.51 i 0.07 

-0.63 i 0.05 

River discharge, Cj) (cm/m 2 s -1 ) 

0.0006 i 0.0001 

— 

0.004 i 0.0004 

— 

.Tangential wind stress, Cy fc (cm/m 2 s“ 2 ) 

— 

— 

— 

-0.04 ± 0.06 

Normal wind stress, (V» (cm/m 2 s“ 2 ) 

"n 

— 

— 


0.85 i 0.07 

St. deviation original record S(ti> (cm) 

7.26 

6.53 

9. 10 

7.80 

St. deviation of residuals R(t^) (cm) 

5.96 

3.44 

4.19 

3.31 

% reduction in variance 

33Z 

72 Z 

79% 



82 Z 


Table 3. Linear regression analysis of selected tide gauge records 
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-latum Bias 

The datum bias terms (C^) merely reflect the 
arbitrary origin of the particular tide gauge 
scale . 

Linear Trend 

This quantity (Cl) may be expected to combine 
secular vertical crustal movements and eustatic 
sea level variations. According to the standard 
deviations obtained, these results are quite well 
determined. Generally, there appears to be a 
high degree of spatial coherence, which is furth- 
er supported by the results for the four omitted 
gauges, despite their lower accuracy. Only Pt. 
Sapin and Pictou were exceptional; the latter 
displaying a negative trend. Overall, the linear 
trend appears to increase fairly smoothly for 
locations towards the south and east. Presumably 
this reflects a regional crustal movement phenom- 
enon, as the eustatic effect should be similar at 
all locations. 

Tides 

Table 3 includes results for both the ampli- 
tude (Aj) and the ascending node date (Nj) for 
the five tidal components. Ascending node date 
is derived from the phase angle 0j of the tide 
by the relationship: 

Nj = B-j + (0j - j ± 2irn) / oij , (5) 

where Bj is the date of the time origin used for 
the particular data series (usually the middle of 
the month immediately preceding the first month 
of data) . The dates given in table 3 have been 
adjusted as necessary by an integral number of 
periods (n) for purposes of comparison. 

Semi-annual Tide . The phase of this component 
is very well determined and spatially consistent. 
A systematic difference in amplitude is apparent 
between those locations affected by river dis- 
charge and those which are not. It seems likely 
that the greater amplitudes found for the gauges 
where river discharge is significant may result 
from poor modelling of that effect. At St. John, 
moderate correlation was found between the phase 
of the semi-annual tide and the river discharge. 

Annual Tide . Amplitudes of the annual tides 
are not as well determined and again the spatial 
coherence seems to depend on the presence, or 
absence of river discharge. There is, however, 
also a strong interaction with the temperature 
effect. For instance, at Halifax fourier trans- 
formations of both the air and observed water 
surface temperature series yielded amplitude 
coefficients for the annual harmonic at least an 
order of magnitude larger than those of any other 
period. High correlation coefficients between 
both amplitude and phase of the annual tide and 
the temperature effect (0.60 and 0.87) respect- 
ively at Halifax) are thus not unexpected. This 
evidence suggests that these two effects are not 
being adequately isolated and, indeed, the phil- 
osophy of attempting to do so might be question- 


able. Ultimately, both the annual tide and the 
strong annual component of temperature variation 
are merely different manifestations of the same 
motivating energy source: the earth's orbital 
motion. Poor coherence in the phase of the 
annual tide may also result from the influence of 
the temperature effect. 

Chandlerian (Polar) Tide . Good spatial con- 
sistency and satisfactory precision are achieved 
here. The amplitudes and phases are also quite 
comparable with results obtained for several east 
coast ports in the U.S.A. [Currie, 1975; Vanicek, 
1978]. 

Lunar Perigee and Nodal Tides . The precision 
attained in evaluating these longer period tides 
is consistent with the total length of the time 
series analysed. Despite fairly weak determin- 
ations, the amplitudes agree reasonably with the 
empirical knowledge of these effects, and also 
with the estimated equilibrium tides [Rossiter, 
1966].. 

Temperature Effect 

Generally, the temperature coefficient (Gp) is 
so weakly determined as to be meaningless. 

Indeed, the negative results obtained at Pt. au 
Phre and Charlottetown would be contrary to 
expectation based on an understanding of the 
physical process involved. At all locations 
except Halifax, air temperatures were used 
directly in the regression model. However, a 
test implementation at Halifax of the water sur- 
face temperature sub-model described above, does 
not noticeably improve the result. In view of 
this outcome, further investigation of the sig- 
nificance of the Interaction between the annual 
tide component and the temperature effects is 
warranted. It seems possible that the lack of 
constraints on the amplitude and phase of the 
annual tide allows undue absorb tion of the data 
variance at that frequency — thus reducing the 
temperature effect — and simultaneously falsify- 
ing the determination of the tidal coefficients. 

Air Pressure Effect 

These coefficients (Cp) are well determined 
and are substantially in agreement with the 
theoretical model and empirical studies of this 
effect. Taking into account the maximum observed 
range of air pressure at the gauge locations, the 
coefficients obtained represent a minimum equiv- 
alent sea level variation of about 10 cm at St. 
John and a maximum of about 20 cm at Charlotte- 
town. 

River Discharge Effect 

River discharge effects are expected to be 
highly variable, depending on the local estuary 
configuration. The coefficients (Cp) obtained 
are equivalent to maximum sea level variations of 
approximately 12 cm and 8 cm at Pt. au Pfere and 
St. John respectively. 
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Wind Stress 


Replacement of the geos trophic wind sub-model 
(eq. 3) by the observed wind sub-model (eq. 4) at 
Halifax appears to have been successful. The 
tangential wind stress effect is effectively 
negligible as anticipated. In contrast, the 
normal component is quite substantial, amounting 
to an equivalent maximum sea level variation of 
approximately 20 cm. Further study of this 
effect will be directed towards establishing an 
empirical relationship between the geostrophic 
and observed wind models, with the aim of improv- 
ing the geostrophic model for use where suitable 
observed wind vectors are not routinely available. 

Residuals 

A worthwhile reduction in variance of the sea 
level record is evident when the standard deviat- 
ion of the residuals R(t^) is compared with that 
of. the original data (table 3, last three lines). 
The same result is depicted graphically in figure 
2, where the plot of residuals may be compared 
visually with the original record. To facilitate 
this visual comparison, the residuals displayed 
in figure 2 (but not the results in table 3) have 
been subjected to a low-pass filter to remove the 
high frequency components with a period shorter 
than about 4 months. 

Spectra of Residuals . Figure 3 illustrates 
the small amount of power remaining in the 
spectra of the unfiltered residuals. If reason- 
able limits are imposed on the band-width, so as 
to be consistent with the data span, there are no 
spectral peaks which exceed about 8% of the var- 
iance. The noticeable remaining variance in the 
Pt. au P&re residuals may relate to the disjoint- 
ed data series or it may reflect the additional 
complexity of the St. Lawrence River environment. 

Conclusions 

While the results presented here must be re- 
garded as preliminary — and any dependent conclus- 
ions consequently tentative — the following sum- 
mary will serve to highlight the strengths and 
weaknesses of the method and hence usefully 
direct further investigations towards realization 
of its full potential. 

(1) A linear sub-model of the air/water heat 
transfer process, incorporating time lag, can be 
constructed empirically with acceptable reliabil- 
ity. However, it would appear that the use of 
water surface temperatures — either derived from 
such a model or actually observed — in the overall 
model for consequent sea level variations, does 
not produce a reliable or coherent thermal coeff- 
icient and does not significantly reduce the var- 
iance of the sea level record. It is suggested 
that a likely cause of this result may be the 
interaction of the highly correlated annual tide 
effect. Even though the coefficient for the 
thermal effect might actually be small — for 
example, Thompson [pers. comm, reported in 
Anderson, 1978] has obtained results which 
suggest that this is the case in the waters near 


U.K. — its weak determination and poor coherence 
indicate a remaining problem in the analysis. 

(2) Use of properly derived observed wind 
vectors produces acceptable wind stress coeffic- 
ients and contributes significantly to the reduc- 
tion of sea level variance. The utility of a 
geostrophic wind sub-model remains uncertain, 
but development of an empirical relationship 
between such a model and the observed wind field 
would seem to be worthwhile. 

(3) A very satisfactory determination of the 
linear trend coefficient is possible and this 
result appears to provide a reliable estimate of 
regional vertical crustal movement. This is 
especially significant in the context of tide 
gauge applications for continental levelling 
datums. 

(4) While the signatures of the polar and 
longer period tides seem to have been adequate- 
ly identified and separated from the sea level 
record, there are remaining difficulties with the 
annual and semi-annual tides. The treatment as 
free parameters of both amplitude and phase of 
these tides, as well as the coefficients of 
thermal and river discharge effects, apparently 
fails to adequately isolate their contributions 
to the sea level variation. Further study of 
this problem is warranted. 

(5) Since modelling of wind stress and baro- 
metric pressure effects is reasonably successful, 
further refinement of the model, by inclusion of 
steric level and ocean current terms for instance, 
might usefully be pursued. This may even assist 
in better resolution of the tidal terms. 

(6) Even in its present preliminary status, 
the modelling process provides sufficient reduct- 
ion of variance in the sea level record for its 
utilization in identification of episodic move- 
ments at the 10 cm level. This facility may be 
further enhanced by investigating the differences 
in the residual variations between pairs of tide 
gauge locations. 
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